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Radar 

Propagation 
Scattering 
Inhomogeneous  Media 

AUSfRACT  M'iMllMtji'  <in  ^Itfr  If  iiri  «»iil  lifrnlilv  hv  l»fi«rA 

The  propagation  of  radar  through  a striated  plasma  has  been  analyzed  on  the 
basis  of  random  walk  photon  scattering,  geometric  optics  refraction  at  a thin 
phase  screen,  plane  wave  scattering  from  a thin  screen  In  the  Fraunhofer  limit 
and  the  observables  to  a monopulse  system  using  a wave  front  integration  Imple- 
mentation of  Huygens'  principle  over  a thin  phase  screen.  The  latter  two  were 
carried  out  for  screens  of  both  large  and  small  phase  variance.  The  scattering 
was  related  to  the  properties  of  the  plasma  for  a variety  of  distribution 
functions  of  strlatlon  sizes.  For  reasonable  distributions,  a similar  
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20.  ABSTRACT  (Continued) 

•'effective'^  strlatlon  size  Is  obtained.  This  allows  the  calculation  of  the 
angular  variance  of  the  scattered  power  which  Is  gaussUnly  distributed  for 
most  significant  cases.  In  the  geometric  optics  regime  the  distribution  of 
density  of  allowable  multipath  rays  Is  g|jMs1anly  distributed  and  the  power 
distribution  Is  Identical  for  all  rays.  The  signal  processor  sees  an  unmodified 
but  attenuated  component  and  a random-phased,  scattered  component  that  leads  to 
Rician  amplitude  statistics.  The  results  of  the  various  approaches  are  all 
compatible  and  furnish  a single  unified  description  of  the  perturbed  signal.  ifL- 


PREFACE 


The  work  ilcscrlbod  hero  was  oriclnallx  performed  In  oarly  1974 
for  the  II. S,  Army  SAIMUlUAUO  System  Command,  and  was  directed  toward  the 
description  of  the  behavior  of  a lllir  radar  operat  ln|j  In  a strongly  scatter- 
ing medium.  Any  iarge-scalo  structure  (tens  of  kilometers)  was  handled 
by  gross  refraction  computations  and  only  the  fine-scale  structure  was 
assumed  to  be  handled  by  the  scattering  analysis,  Most  aspects  of  the 
analysis  covered  hero  are  precisely  applicable  only  in  the  limit  that 
plasma  structural  dimensions  are  small  compared  to  a Fresnel  sono.  This 
assumption  Is  identical  to  the  assumption  that  the  receiving  antenna  is  In 
the  far  cone  of  tlic  scattering  medium  when  It  Is  considered  to  bo  a trans- 
mitting source.  Uccent  research  has  shown  that  extreme  care  Is  necessary 
In  the  c.'stcnslon  of  theory  to  higher  frequencies  where  the  structural 
assumptions  may  not  bo  valid. 
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PROBLEM  STATEMENT 


Tlio  ilctoiiiitlon  of  luiclciir  \voa|ions  iiliovc  about  81)  kjii  results  in 
cxtonsLvo  rcjjlons  of  i onl  iiit  ion  ■— jiart  1 al  ly  lonkoil  plasmas.  Dotonatlon 
Liulucoil  winds  iiiiiy  blow  this  plasma  across  the  jtcoma(.',notlc  field  in  a way 
tluit  I’osults  in  install!  1 itles  in  the  lonliod  component,  Such  install i 1 1 1 ies 
produce  striations  e.\tended  alonn  the  Hoomagnetlc  field  lines,  Those  strla- 
tions  have  been  seen  on  all  nuclear  tests  at  very  hluli  altitudes. 

If  radars  attempt  measuremonts  through  striated  renlnns,  the 
wave  front.s  are  liroken  up,  resulting  in  scattering  and  multipath  prnpanatlon, 
The  quantification  of  the  effects  produced  l)y  striations  has  continued  in 
the  effort  reported  here,  Specifically,  additional  analysis  has  unified 
prior  work,  extended  models  of  the  effects  to  weak  scatterln}’  cases,  removed 
some  arbitrary  assumptions,  identified  limits  of  applicability  of  previous 
models  and  considered  tlu’  scnsltivlt>’  to  \'arlous  assumptions  aliout  plasma 
sii'uciure.  t!u treat  thcoreiic.il  pliisma  aatilysis  .ciolds  little  information 
about  striation  ioni?,atlon  profiles  or  ,1  i st  r 1 Init  1 on  of  striatlon  sittes. 
tionsciiueni ly , unless  a relative  l nscitsi 1 1 v 1 t.c  to  those  plasma  parameters  is 
demonst rtited,  modellitK  of  the  derived  effects  will  luive  limited  utillt,v. 

Several  approaches  to  the  scattering;  problem  will  be  considered. 
I.ach  liiis  its  own  region  of  validity  and  limitations  and  is  useful  to  Illus- 
trate particular  aspects  of  the  scattering,  Specifically,  the  first  approach 
will  be  a photon  scattcrlns;  anal\’sls  that  is  valid  for  conditions  under  wlilch 
photons  are  I'luise- 1 ncoherent , Subsequent  iqiproacltes  will  use  the  a|iproxima- 
tion  that  the  1 nhomo>;eneous  medium  is  relatively  localized  and  can  be  aiiproxl 
mated  liy  comblniin;  tl\e  integrated  plasma  Induced  phase  shift  in  a slitule, 
lliln,  phase  slilftiny,  screen,  tleneral  cliai'actei'l  st  1 cs , valid  when  Keoraetrlc 
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optics  approximations  apply,  will  bo  derived.  A physical  optics  analysis 
using  lluygLMis'  principle  will  obtain  a genornlizod  angular  power  spectrum 
for  cases  of  cither  small  or  largo  phase  perturbations,  following  these 
basic  scattering  anaiysos,  the  interpretation  of  the  multipath  signal  by 
a phase  comparison  monopulse  system  will  be  developed,  the  sensitivity  to 
alternate  models  explored  and  scattering  In  tei'ins  of  plasma  parameters 
summarl zed . 


THEORETICAL  AND  EXPERIMENTAL  BACKGROUND 


2,1  Plasma  Induced  Phase  Shift 

for  radar  waves  (which  have  frcquenc Ic.s  much  greater  than  cither 
electron  collision  or  gyro-fretiuenclos) , the  Index  of  refraction  of  a plasma 
Is  given  by 


u*  ■>  1 


where  go  ' the  permittivity  of  free  space,  n^^  is  the  electron  concen- 
tration, e and  m are  the  electronic  charge  and  mass  and  uj  Is  the 
radar  angular  frequency.  I'enote  a critical  electron  content 


so  that 


n^  = u’  mco/c*, 


U*  " I - n^yn^. 


lk)r  conditions  of  Interest  to  radar  propagation,  usually 
n^^  « n^.  tionsciiuontly , 
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This  apiu'oximat Uin  will  he  used  In  the  reimiiiulcv  of  this  report. 


The  electric  vector  of  a plane  wave  propapatliiR  in  a plasma  is 

Riven  by 


(x.t) 


= l-j  0 


(S) 


Tills  iiieiuis  the  phase  Is 


ij'  ) ■-  wi  - (o.\p/c 


i<0 


Compared  to  a wave  in  the  absence  of  plasma  01“  l.)»  fh*-'  I'hase  advance  per 
unit  path  Is 


..  TiLllJil  ^ 11'.  .S'.. 

d.\  cu  c in 


(7) 


2.2 


Str1at1on  Profile 


f 


ior  most  of  tile  detailed  analysis  in  this  report,  tlie  individual 
strlatioiis  will  lie  assumed  to  he  c>l  liulr  lea  I !>■  symmetric  and  have  an  electron 
concentration  iiroflle  that  Is  Rausslan.  There  Is  no  e.xper Imenta I data  to 
support  or  contradict  this  assumption.  If  the  profile  were  controlled  e.\- 
cluslvely  by  diffusion  and  the  diffusion  coefficient  were  constant,  It 
would  he  Ruussian.  However,  the  real  rationale  for  usiiiR  this  model  is 
that  it  renders  the  mathematics  tractable.  The  justification  for  accep- 
tance of  the  results  is  that  where  alternate  assumptions  can  he  used,  the 
results  are  not  sIriiI f leant  ly  altered.  Some  of  these  cases  will  he  pointed 
out  in  the  follow! hr  discussions. 


'I’lio  olc'ftron  concentration  will  bo  Riven  by 


whore  no  is  tlie  iixlal  electron  coneentrnt  ion , r Is  the  i.llstiince  (-‘rom 
the  axis  aiul  o is  the  characteristic  dimension  of  the  strlatlon,  lUstri- 
luitlons  of  sl;es  will  be  considered  and  no  i«ay  be  a function  of  o. 

2.3  Experimental  Strlatlon  Size  Distribution 


The  only  expor imenta  I data  oi\  strlatlon  sl:e  distributions  for 
miclear  events  comes  from  photographic  data.  Thus,  It  is  limited  to  regions 
of  optical  emission,  the  highly  disturbed  or  fireball  regions.  There  is 
theoretical  basis  for  believing  that  strl.atlons  .also  form  outside  high  alti- 
tude fireballs.  It  must  be  remembered  then  that  the  experimental  datii  Is 
obtained  from  onl.v  one  type  of  unstable  region  and  Is  not  necessarily  rep- 
resental ive  of  all  such  regions. 

Pertinent  photograpliic  data  have  been  analysed  by  W.  C'hesnut^'\ 
SRI,  and  the  data  can  be  reasonably  well  fit  by  gaussi.-in  strlations  with 
constant  oi\-axls  concentration  and  a site  distribution 


('o  ) 


Some  of  the  data  were  best  fit  by  the  sujH'rptis  1 1 Ion  of  two  distributions, 
each  of  the  form  of  lUpiat ion  P.  Also,  the  quality  of  the  data  is  insufficient 
to  define  the  exponent  of  the  oq/o  term  very  well — f,  or  might  be  as  good 
as  ii.  but,  use  of  the  sixtii  power  iiuiKes  some  desired  integrals  closed,  well 
ki\own  ana  hole  functions. 
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2.4  Strlatlons  Outside  Fireballs 

Theoretical  analysls^^^  has  shown  that  where  the  direction  of  the 
explosion  induced  wind  is  both  generally  crossing  the  geomegnetic  field 
lines  and  in  the  direction  of  decreasing  electron  content,  the  ionized 
fluid  flow  is  unstable  and  should  develop  perturbations  (and  possibly 
strlatlons).  Spatial  structure  of  characteristic  wavelength  greater  than 
a minimum  cutoff  wavelength  should  grow. 

If  the  environment  is  assumed  to  be  composed  of  the  superposition 
of  a large  number  of  such  Fourier  components,  the  spatial  amplitude  of  which 
is  limited  to  the  spatial  wavelength,  the  number  of  strlatlons  tends  to  fall 
off  as  an  inverse  power  of  the  striation  size.  Below  a minimum  size,  the 

population  should  be  nearly  zero.  The  amplitude  of  the  disturbance  may  be 

f 31 

a function  of  the  striation  size.  Previous  HARC  analysis  assumed  that 
the  concentration  perturbation  n^  was  proportional  to  the  striation  wave* 

length  (due  to  interchanging  air  over  a region  a wavelength  long)  and  that 

♦ 

the  appropriate  power  law  was  the  inverse  fourth. 

It  is  necessary  to  relate  cutoff  wavelength  and  minimum  striation 
size.  In  the  HARC  work,  the  sine  wave,  Fourier  components  were  approximated 


The  number  of  strlatlons  of  a given  size  per  unit  area  perpendicular 
to  the  magnetic  field  was  taken  to  be  Inversely  proportional  to  the 
size  and  proportional  to  the  allowed  mode  density  on  a circular  ring. 
The  mode  density  is  then 


dn  2tiR 

3x  “ ir  • 


JX  IK 


by  tt  series  of  uniformly  spaced  gaussian  striations.  A reasonable  approxi- 
mation was  obtained  if  o*  » X*/22  . Alternatively,  if  the  relationship  is 
based  upon  either  the  sine  wave  or  a series  of  gaussian  striations  which 
have  equal  inuxiinuin  concentrations  and  produce  the  same  rms  scattering 
coefficient,  the  relationship  should  be  o*  ■ ■ XVl7.B  , Conse- 

quontly,  u ratio  of  the  order  of  20  is  reasonable. 


PROPAGATION  ANALYSIS 


3.1  Random  Walk  Approach  to  Scattering 

Inasmuch  as  the  density  of  striations  is  relatively  groat,  a 
particular  ray  will  interact  with  several  and  exhibit  motion  that  cun  be 
treated  as  u random  walk  (in  angular  deviation).  The  bonding  of  a ray  in 
an  inhomogeneous  medium  is  given  by 


where  ^ indicates  u spatial  derivative  perpendicular  to  tho  direction  of 
propagation,  (Because  tho  striations  are  extended  along  tho  geomagnetic 
field,  the  problem  can  be  treated  us  two  dimensional  — all  bonding  being  in 
tho  plane  normal  to  tho  goomugnotic  field.) 

liquations  7 and  8 can  bo  used  to  obtain  tho  deviation  duo  to  one 
gaussian  striation.  If  tho  bending  is  groat,  system  operation  will  fail. 

It  follows  that  for  coses  of  interest,  the  ray  path  will  be  nearly  u straight 
lino  and  tho  deflection  can  bo  approximated  by  integrating  the  incremental 
bonding  along  a straight  line.  Lot  a particular  sight  lino  poss  a distance 
b from  the  striation  axis.  Therefore, 


6ecb)  " nX,  ^al  • 


The  integration  and  ^ are  independent,  and  their  order  can  be  reversed 


“ H A / 


art 

C 


at 

r 

J 2ti. 


b®  + (x  - ) * sin*  a 


whoro  the  jiolnt  of  closest  approach  of  the  strliition  axis  and  the  sight  path 
is  at  and  a Is  the  angle  between  the  sight  path  and  the  strlation 
axis,  The  terminals  of  the  path  are  sufficiently  far  from  the  strlation 
axis  that  there  Is  little  error  produced  by  placing  the  limits  of  Integra- 
tion at  1^’^. 

siuni  ■ -bV2o* 

‘ T,rjrni^  iir 

c 


In  this  form  the  perpendicular  derivative  Is  the  derivative  with 
respect  to  b or 

n.VSiT  . -b*/2ci* 

W(b)  - ^ e . (16) 

After  a random  walk,  the  finol  value  tends  toward  a gaussian  distri- 
bution with  a variance  equal  to  the  sum  of  the  squares  of  the  individual 
steps.  The  value  of  the  mean  of  this  sum  can  be  obtained  by  multiplying  the 
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square  of  tlie  incremental  deviation  (Equation  16)  by  the  probability  of 
having  n atriution  of  site  o at  a displacement  b and  integrating  over 
all  0 and  b.  Thus, 


vr 


II 


2n"  sin^'  ct 
c 


I 


b ■ 


b* 


-b®/c5® 

e l’(n)  sin  a mb  do  db 

(17) 


where  m is  the  concentrntion  of  strlations  on  a plane  normal  to  their  axes, 
L Is  the  length  of  the  path  through  tho  striated  region  and  P(o)  is  the 
probability  density  function  of  o. 


0 


2 

0 


« 


1 n^,  sin  a J 
0 • 0 


(18) 


We  will  use  tho  notation 


“T"? 

"n  " 


(It 

I 


n)  o'Tio)  do  . 


0 • 0 

Thus,  Equation  18  can  bo  written 


(U>) 


0 


2 

0 


4 n*  sin  a 

L' 


(20) 


This  equation  allows  tho  computation  of  tho  variance  of  tho  scattering,  and 
after  numerous  scattering  events,  the  distribution  of  power  tends  to  a 
gausslun  distribution.  The  guussian  character  is  the  result  of  the  multi- 
plicity of  events,  not  the  distribution  of  angles  due  to  a single  interaction 
Consequently,  it  is  not  a function  of  tho  electron  profile  in  tho  striotlon. 
In  fact,  the  scattering  from  a gausslun  rod  is  far  from  gaussian  but  exhibits 
discrete  maximum  scattering  angles  each  side  of  zero  and  has  maximum  proba- 
bility at  tho  extremes  of  tho  distribution. 
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The  restriction  on  the  applicability  of  this  analysis  is  that  all 
of  the  potential  scatter  paths  contributing  energy  in  a certain  direction 
must  be  phase  incoherent  in  order  that  interference  patterns  do  not  exist. 

As  will  be  discussed  later,  for  rms  scattering  in  excess  of  a few  milll- 
radians,  for  typical  structural  dimensions  and  UHP  radars,  phase  Independence 
should  be  assured, 


In  one  sense,  the  random  walk  analysis  is  the  most  general  treat- 
ment in  the  strong  scattering  case.  It  is  applicable  even  for  environments 
that  produce  crossing  of  rays  within  the  medium. 

3.2  Phase  Screen  Approximation  to  Strong  Scattering 

The  phase  screen  approximation  to  refraction  and  scattering 
assumes  the  scattering  medium  is  sufficiently  localized  that  its  effect 
can  be  approximated  by  considering  that  all  the  phase  shift  experienced 
by  a ray  traversing  the  medium  can  be  attributed  to  a thin  screen  in  the 
middle  of  the  medium,  It  also  implies  that  the  amplitude  is  constant  over 
the  screen,  This  section  will  cover  parameters  describing  the  screen  and 
the  relationship  to  volume  characteristics  of  the  medium. 


Usually  the  characteristic  assigned  to  the  screen  is  the  integral 
of  the  phase  shift  on  a line  parollel  to  the  undisturbed  ray.  That  is,  if 
the  radar  is  at  - Rp,  0 and  the  target  is  at  R,j,,  0, 


. A d^(x,y) 

J ox 


In  the  limit  that  the  phase  is  the  result  of  a random  process  and 
irregularities  fill  Fresnel  zones,  the  angular  deviation  Is 


L.l, 


Close  to  the  phase  screen i these  conditions  are  valid  because  Fresnel  zones 
are  small  compared  to  the  typical  strlatlon  sizes  of  many  hundreds  of  meters. 
Under  those  conditions,  It  is  also  possible  to  determine  the  power  in  the 
multipath  rays  by  considering  the  region  over  which  a wave  Is  scattered 
coherently,  l.e.,  the  extent  over  which  the  "scattered  ray"  varies  by  only 
X/2.  This  is  controlled  by  the  second  derivative  of  the  phase  shift. 

Lot  us  first  develop  expressions  relating  the  phase  shift  to 
plasma  parameters  and  then  determine  the  scattering.  Then  the  scattering 
will  bo  discussed  in  terms  of  the  characteristics  of  the  energy  at  the 
radar. 

The  mean  integrated  electron  content  associated  with  the  striations 
is  obtained  by  forming  the  product  of  the  probability  that  a striation  axis 
lies  at  some  Impact  parameter  and  the  integral  of  electron  content  for  that 

striation  and  impact  porameter,  and  integrating  the  product  over  all  Impact 

* 

parameters. 

« « « * x V siV  g) 

^■///  rtj  e slnamP(a3  L dx  db  do  (23) 
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The  contribution  of  the  striations  to  the  local  electron  concentration  is 
just  N/L.  The  striations  normally  are  superimposed  upon  a constant  (or 


The  limits  of  on  the  x integral  reflect  the  assumption  that 
the  striation  sizes  o are  small  compared  to  L so  that  little 
error  results  in  extending  the  limits  of  integration  rather  than 
truncating  them  at  their  actual  value. 


slowly  varying^  eoncontration  no  that  may  be  considerably  greater  than 
the  mean  strintlon  contribution.  The  total  mean  electron  concentration  is 

» nij  + N/L  . (25] 

A parameter  some  plasma  physicists  feel  they  can  compute  is  the 
local  mean-square  electron  concentration  fluctuation  due  to  striation 
structure. 

2 ~ — * 


ww 


? 0 ® 2nr  mt’(CJ)  dr  do  - ({T/L)* 


■ flin  n8  0*  - (N/D*  (27) 

In  the  remainder  of  this  report  the  term  (?3’/L)*  will  be  dropped,  an  approxi- 
mation valid  for  striations  sufficiently  separated  that  the  primary  contri- 
bution to  0*  is  from  regions  without  striation  overlap  (i.o.  4'?rm"o*  1)  . 

The  mean  phase  shift  through  the  screen  is  unimportant  to  the 
problem  being  addressed  but  could  be  easily  obtained  from  Expression  25. 

In  the  following,  the  direction  of  primary  pi'opagatlon  (the  sight 
path)  la  nearly  parallel  to  the  x-axis,  the  striation  axes  lie  parallel  to 
a line  in  the  x-z  plane  that  makes  an  angle  a relative  to  the  x-axis, 
and  the  coordinate  y is  perpendicular  to  both  the  sight  path  and  the 
striation  axes.  Moreover,  any  fixed  background  of  electron  content  not 
associated  with  striations  has  been  suppressed  because  it  does  not  affect 
the  scattering  and  is  accounted  for  in  refraction  calculations. 


The  variance  of  phase  ehift  is  important  because  it  determines 
the  magnitude  of  phase  shift  fluctuations.  Consider  a particular  striation 
with  its  axis  at  y ■ h^.  The  integrated  phase  shift  along  the  line  y ■ 0 
is 


6(|i^ 


m 


J —2^1 


dx 


(28) 


yTfT  ir  Oj  -b*/2o| 
" V sin  a ® 

\/Tn  tr 

- xi^rnn;  >’i  ' "i  “i- 

The  phase  shift  variance  is 


(29) 


aj  . - ?*  (30) 

whore  <t)  is  u sum  over  all  strlntions,  or 

oj  ■ C*  [ ^ P^  exp  (-b|/2o|)  j*-  . (30 

The  bar  designates  an  average  over  all  possible  configurations  of  atria- 
tioiis.  This  con  bo  rewritten 

'i 

oj  ■ p1  oxp(-b*/o|)]  j 

' :) 

+ C’  ^ L Pj  exp  ( -bj/ZoJ  > b*/2o*  ) - T*  . (32)  j I 

The  assumption  that  the  striations  are  randomly  located  implies  that  h^  | f 

is  independent  of  bj  so  that  the  mean  of  a product  of  functions  of  i ; | 

and  j osiuals  the  product  of  the  individual  means,  or  i j 
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21  Pi  «xp 


♦ C»  22  Pi  «XP  (-bi/2(j5)E^  Pj  «xp  ( -b*/20»)  C33) 

The  last  term  in  Equation  33  is  made  up  of  two  factors,  the  first  la  7 
and  the  second  ? (1  - ^ ) where  k Is  the  number  of  striatlons.  If  k 
is  larg(  . th > 1/k  will  be  negligible, and  it  will  be  so  assumed  here. 

Thus 


0*  - C»  L pj  oxp  ( .b*/oJ  ) • 


(34) 


For  large  ensembles,  the  summation  and  averaging  can  bo  approxi- 

•it 

mated  by  integrating  over  the  spatial  and  size  distribution  functions 


0 


2 


00  00 

c;/[ni  c 


oxp(-bVa*)  ml  sin  a PCo)  db  do 


(35) 


■ 


ml 

n^  sin  d 


n?  o’. 


(36) 


The  factor  sin  a occurs  In  Equation  35  because  the  concentration  of  strl- 
ations  In  the  xy  plane  is  m sin  a. 


j 

-a 


* The  result,  Equation  36,  can  be  combined  with  Equation  27  to  give 

a*  ■ — -p --y  where  o » nJoVn?  o“,  This  is  the  same 

™ ‘‘c 

form  derived  by  Bruniloy^^^  for  a thick  scattering  region. 
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Another  Function,  useful  because  It  determines  the  scattering, 
is  the  autocorrelation  function  of  the  phase  shift.  It  can  be  computed 
In  a manner  similar  to  that  used  for  the  variance  of  phase.  The  nor- 
malized autocorrelation  function  Is  defined  as 


R (O 


<t>(y'j  - ? JL»(y*C)  - 


^(y)  jiCyH)  - 


(37) 

(38) 


(39) 

( R will  be  used  to  denote  various  autocorrelation  functions  but  when 
unsubscrlptod  will  denote  the  phase  autocorrelation  function.)  By  a sequence 
of  steps  analogous  to  those  in  Equations  31  through  34,  Equation  39  can 
be  reduced  to 


Again,  the  average  of  the  sum  can  be  approximated  by  an  integral  over  space 
and  strlotion  parameter. 
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I 


IB 


i 

i 


R(5)  . , .,..p5,J"ii..^.  f o’  cxp(.  -ii)  P(a)  do 

X’n*  3ina  o*  *'*  \ 4o*/ 

c (p 

Oi',  if  the  value  of  a?  is  substituted  from  P.xijression  36, 


(413 
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1U53  - 


n;  O' 


(42) 


How  does  this  relate  to  ray  deviation?  Consider  the  general 
normalized  autocorrelation  function  whore  for  simplification  In  writing 
it)  is  taken  as  the  deviation  from  the  mean. 


11(5)  . Urn  j r^Uo’)  ♦(y  + C)  Jy!  /2yiO» 

yiH«(-'.yj  ) 


I /2v.  a* 


(43) 


The  second  derivative  of  11(5)  ut  5"d  is 


II"  (OJ 


lim  I i|)(y)4)"  (y)  dy  | /2y^a^ 

( •'-xi  ) 


(44) 


Integrate  this  by  parts  once  to  obtain 


U"(0)«  11m 
'i 


y,-*^  I 2y,  o’ 


^ r<^'(y}l*  dyj.  (45) 

) \ 


If  <))  and  4^'  remain  small  over  the  range  of  interest,  the  first  term  of 
the  right  hand  member  becomes  small  compared  to  the  second.  Moreover,  the 
second  term  is  the  ratio  of  the  variances  of  the  derivative  of  phase  and 
the  phase.  Consequently, 
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(46) 


oj.  --°JK"C0)  . 


( ^'  is  nssumod  to  have  zero  mean). 


The  deviation  is  directly  rointed  to  the  dorivative  of  phase 
shift  CEiiuiition  22).  Thus, 


X^o 


4ir' 


411* 


R"(0)  . 


(47) 


As  used  hero  is  the  variance  of  ray  bcndini}  or  seuttorinR  nnRlc.  In 
later  analyses,  It  loses  that  physical  significance  hut  is  still  a useful 
characterization  puranictoi’  for  the  plasmas.  This  cun  be  used  with  I'.xprossion 
41  to  give 
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,3/! 
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ml. 


® 4 sin  u n?, 


ni  0 


(48) 


which  is  idonticai  to  the  result  given  in  Ikjuation  20. 
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Ray  Optics  and  Stationary  Phase  Analysis 
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3.3.1  Scattering  as  Viewed  by  the  Radar 

The  iieonietry  of  the  thin  phnso  scroon-radar-target  in  shown 
In  I’lgure  1.  Up  until  now  the  scattering  angle  0 has  boon  discussed, 
tt  can  be  scon  that  the  angle  of  scatter  seen  by  the  radar  ip  is 

4/  ■ 9 R.j,/(r+R.j.)  (49) 

if  the  angles  are  .small. 

In  the  strong  scattering  case  (tho  phase  scroon  having  many 
radians  phase  shift),  unless  there  is  n region  of  the  screen  for  which 
the  scatter  path  phase  shift  is  constant,  the  scattered  energy  will 
add  incohorontly  ond  contribute  a nolse*like,  random  phase  background. 

Tho  scattered  waves  from  regions  of  constant  or  stationary  phase  will 
coiiiblno  coherently  and  form  discrete  scatter  ray  paths. 


Figure  1.  Scatter  geometry. 
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The  total  phase  advance  Including  geometric  and  phase  screen 
contributions  for  a path  penetrating  the  screen  at  y is  (for  small 
angles,) 


whore 


+ <t>  . 


U'*  ■ r 


Fermat's  principle  leads  to  the  condition  that  rays  between  two  points 
(radar  and  target)  propagate  over  paths  of  local  extreme  phase  deviation, 
i.e,,  where  the  derivative  of  phase  is  zero.  These  are  culled  points  of 
stationary  phase.  If  wo  set  the  derivative  of  oqunl  to  zero,  we  get 

T'  D'l ' ■ . (52) 

where  f has  been  used  to  represent  AH/2'it,  Hquutions  40  and  52  represent 
the  same  relationship  if  9 “ ^ < which  wos  previously  derived  us  the 

scattering  angle  at  the  phase  screen. 

The  intensity  from  one  of  these  regions  is  determined  by  the 
extent  of  the  coherency  region.  The  ampiitudo  varies  ns  the  linear  extent, 
the  power  us  the  square  of  the  extent.  (Quantitatively,  the  extent  will  be 
about  that  for  which  j Ati...  1 . That  is. 


I y ^ A 2 

<t>.].Cy)  - <ti.[.(y^)l*  n ■ - — Ay  + (|)"  iVj)  , 


where  the  phase  shift  has  been  expanded  os  a Taylor  series  truncated  after 
the  quadratic  term  and  Ay  ■ y - y^.  Algebraic  manipulation  of  this  using 
the  value  of  <[)'  from  Equation  52  produces 


Ay*  - 2Tf  -I . r (yi)| 


For  no  atriations  or  refraction,  if"  Is  zero  everywhere.  The 
use  of  this  normalizes  Equation  S4  to  the  undisturbed  power  gain. 


JL  . 1 

I'o  1 - fif"  (y^) 


Although  Equation  55  indicates  that  P/Po  could  go  to  infinity 
if  fif"  « 1 , this  is  u fallacy  introduced  by  dropping  high-order  deriva- 
tives from  Eiiuation  55.  If  tho  structure  of  the  phase  screen  has  a charac- 
teristic dimension  Si  and  scattering  is  strong,  it  is  unlikely  that  the 
phase  is  "stationary"  for  more  than  a distance  H/IO  , If  this  is  used  to 
estimate  an  ujqior  limit  on  the  power  gain, 


2boW 


YoTfXR 


The  analysis  in  this  section  does  not  require  that  the  regions  of 
stationary  phase  be  as  large  as  a Frer-nel  zone  (Jlj.  ■ Ml),  but  it  is  apparent 
from  Equation  56  that  the  power  gain  on  a ray  path  is  the  square  of  the 


extent  of  coherence  measured  in  terms  of  the  first  Fresnel  zone  size.  Con- 


sequently, if  the  structure  is  fine,  most  of  tho  energy  is  either  scattered 
out  of  tho  beam  or  combines  noncoherent ly. 


Distribution  In  Phase,  Angle  and  Power 


The  random  walk  analysis  leads  one  to  believe  that  the  angular 
scattering  distribution  sliould  tend  toward  a gaussian  distribution.  Does 
tho  .stationary  phase  analysis  give  tho  same  answer?  And  if  so,  can  one  draw 
conclusions  about  tho  distribution  of  phase  shifts  and  amplitudes?  The 
anuvor  is  yes. 


If  a straight  path  through  the  medium  intersects  many  randomly 
placed  striations,  the  number  Intersected  will  have  a Poisson  distribution, 
which  approximates  a gaussian  for  large  numbers  of  intersections,  Con- 
sequently, in  this  limit  ifCy)  is  gausslanly  di.stributed. 


Call  its  normalized  autocorrelation  function  R(C)  • This  means  that  if 
the  phase  screen  is  randomly  siunpled  at  points  sufficiently  spaced  that 
RCC}«1  I the  values  so  measured  will  form  a s<^usslan  distribution  with 
variance  , as  already  computed. 

Tor  such  a random  gaussiun  sequence,  samples  spaced  Ly  apart, 
satisfy  the  relationship 

t^(y•*•Ay)  ■ <|i(y)  R{Ay)  + A-WHKy)  , (57) 


whore  g is  a guussiun  variate  from  u unit  variance  distribution,  This 
allows  the  computation  of  the  derivative  of  iji. 


Ay*0 


(SB) 


To  evaluate  this  derivative  requires  knowing  IHO  only  near 
the  origin.  As  wins  proviou.sly  shown  (liquation  46),  for  well  behaved 
sequences,  the  autocorrelation  function  Is  parabolic  near  the  origin. 
Thus , 


«({;)  - 1 


a5'/2  , 


(59) 


a Just  being  used  for  for  convenience. 


The  combination  of  liquations  57,  58  and  S9  yields 


r - lim  + *^>'2  BCT^  ■ 

Ay-*’!)  5y 

- *ZS-  go , . o.'g  . 


(60) 

(61) 
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It  Is  evident  that  the  gausslan  character  of  the  0 distribution 
is  adociuate  to  guarantee  a gausslan  distribution  of  41',  Intuitively 
we  suspect  that  the  conditions  necessary  to  produce  a gausslan  4) 
distribution  are  those  necessary  to  produce  a gausslan  distribution 
in  the  random  walk. 

The  other,  very  significant  implication  of  liquation  61  is 
that  the  slope  of  the  curve  4>'  is  independent  of  the  value  of  41  , 
This  means  that  at  points  of  stationary  phase,  the  phase  shift  0 is 
not  related  to  tho  deflection  which  is  determined  by  0'  but  rather  is 
independently  guussianly  distributed  with  its  own  variance  . Con- 
sequently, so  long  us  0^  ii*  1 , the  phases  of  the  various  rays  are 
indeed  uncorrelutod. 

The  proper  gausslan  statistics  for  0 and  correlation 
between  siunplos  are  maintained  for  three  equally  spaced  samples  if  tho 
second  is  determined  by  Equation  57  and  the  third  is  determined  by 

. . 2H‘(Ay^»(2tyl  ■ ZBMty)  ■ , (52) 

whore  U ■ 1 - H*(Ay)  (Keforonce  5). 

This  can  bo  used  with  tho  autocorrelation  function  to  compute  tho 
second  dorivative 

4,H(y)  . Urn  ; (63) 

Ay>0  ^ 

however,  torina  of  hiuhor  order  than  must  be  carried  In  tho  auto- 
correlation function  In  order  to  oroporlv  represent  the  first  signifi- 
cant term  In  the  rndlcnl  of  liquation  62. 
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Tho  necosaary  expression  for  the  autocorrelation  function 


11(0  - 1 - a5*/2  4 bo/2‘1  . 


The  coefficient  b is 

il'Ml  tv 
0 

This  can  bo  obtained  by  takinu  tho  dorivotive  with  respect  to  £ of 
liquation  41  iiiul  then  setting  5 ■ 0.  Thus 

«J,‘'(0)  . . (66) 

nj  0 

By  taking  the  fourth  derivative  of  tho  general  equation  for  tho  auto- 
correlation function  (liquation  42)  and  integrating  by  parts  twicoi 
we  also  obtain 

k!/(01  . (67) 

Tliese  will  allow  the  determination  of  tho  second  derivative  correlation 
when  liquation  63  is  evaluated. 

Tho  evaluation  of  liquation  63  using  expressions  57,  62  and 
63  is  straightforward  algebra,  but  tedious.  The  result  is 


(()"  ■ . mji  + v1)-a^  go^ 

• • ^ * J'  5^ 
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The  Implication  of  Equation  69  Is  that  the  second  derivative 
of  1^1  which  determines  the  strength  of  a multipath  ray,  la  correlated 
with  the  phase  shift  on  that  ray  and  has  a random  component  as  well. 
Because  there  Is  no  correlation  between  the  phase  shift  and  the  first 
derivative,  it  follows  that  there  Is  also  none  between  the  first  and 
second  derivative  (a  result  of  its  dependence  on  41  Itself).  The  same 
independence  can  be  proven  because  the  first  derivative  Is  a gaussian 
sequence  similar  to  (p  according  to  Equation  62,  and  consequently 
its  derivative  is  uncorrelated  with  itself. 

The  deductions  from  those  correlation  analyses  are  that  at  any 
point  on  the  phase  screen,  1^'  Is  gausslanly  distributed,  and  hence  the 
probability  of  a stationary  phase  ray  occurring  at  a particular  value  of 
y is  a gaussian  function  of  y In  order  that  relationship  Equation  52 
is  satisfied.  The  phase  screen  phase  shift  along  a ray  path  Is  also 
gausslanly  distributed  and  uncorrelated  with  the  scattering  angle.  The 
second  derivative  of  the  phase  shift  which  determines  the  strength  of  the 
signal  propagated  along  a ray  is  also  gausslanly  distributed,  being  the 
weighted  sum  of  two  gaussian  variates  — 4 and  g . The  Intensity  of  a 
ray  is  uncorrolated  with  the  angular  deviation  because  4)'  and  4>" 
are  uncorrolated  — all  possible  paths  have  the  some  energy  probability 
distribution  function.  Consequently,  because  the  number  of  paths  varies 
gausslanly  with  angle,  the  total  power  scattered  varies  gausslanly  with 
angle. 


3.3.3  Density  of  Multipath  Returns 

Previous  analysis  has  shown  that  when  multipath  scattering 
conditions  exist.  If  consideration  Is  taken  of  all  the  combinations 
of  outgoing  and  return  paths,  the  numbor  of  combination  rays  becomes 
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very  large.  They  are  typically  spaced  at  effective  range  separations 
of  a fraction  of  a meter  over  a region  with  a variance 


If  the  radar  is  not  an  extremely  short  pulse  system,  many  signals  that 
appear  randomly  phased  will  interact  In  the  monopulse  processor  as 
was  assumed  In  the  original  IIAKC  analysis.  Consequently,  that  analysis 
which  assumed  the  combination  of  randomly  phased  vectors  In  the  sum  and 
difference  channels  of  u monopulsc  processor  Is  Indeed  valid  for  the 
strong  scattering  case. 

The  sensitivity  of  radar  output  to  variations  in  oiivlron- 
mentnl  models  will  he  discussed  later.  The  extension  of  the  theory 
to  the  generalized  physical  optics  case  will  be  discussed  next. 

3.4  Huygens'  Wave  Theory  and  Angular  Power  Spectrum 

Huygens'  wove  formulation  is  that  at  u particular  Instant  in 
time  one  may  consider  a now  sphoricully  divergent  wave  start.s  at 
every  point  on  a wave  front.  The  superposition  of  these  wavelets  is 
the  solution  to  the  wave  equation  at  n future  point  in  time.  This 
formulation  is,  in  fact,  the  basis  for  the  analysis  of  the  previous 
section  where  the  regions  of  stationary  phase  Imply  that  the  wavelets 
add  coherently. 

In  this  section,  the  Initial  wave  will  be  assumed  to  bo  a 
plane  wave,  there  will  be  a phase  shift  Introduced  by  a thin  phase 
screen,  the  angular  region  of  Interest  will  be  limited  to  nearly 
forward  scatter  so  that  the  wavelet  amplitude  Is  Independent  of  angle, 
and  the  angular  spectrum  will  be  analyzed  far  from  the  screen  so  the 
range  to  the  screen  Is  nearly  Independent  of  position  (from  the 
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standpoint  of  power  divergence,  not  phase  shift).  The  angular  power 
spectrum  will  be  obtained  for  limiting  cases  of  large  and  small  phase 
perturbations  and,  for  one  strlatlon  distribution,  for  all  phase 
perturbations. 

(7) 

Booker^  has  shown  that  the  angular  distribution  of  power 
scattered  by  a piano  wave  Incident  on  a plane  phase  screen  can  be 
obtained  as 
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whore  is  the  autocorrelation  function  of  the  electric  field 

and  Po  is  the  integral  of  P(6)  over  all  6,  equation  71  Is  valid 
In  regions  whore  sinB  » 0 — l.e.  most  energy  scattered  at  small 
angles. 

3.4.1  Autocorrelation  Function  of  Electric  Field 

The  autocorrelation  function  of  the  electric  field  is 


- B(yy  E* (y+t)  , (7: 

where  the  asterisk  denotes  the  complex  conjugate  and  the  field  strengths 
are  normalized  to  unity.  Consequently, 


RljCC)  ■ oxp 


The  phase  shift  at  y-t^  can  be  expressed  in  terms  of  the  phase  auto- 
correlation function  as  given  In  liquation  57. 
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Rp(e)  ■ exp  i i «(y)  [l-HK)]  - 1 v^l-RCO*  g j . (74) 


Inasmuch  as  ({)(/)  and  g are  independent,  the  average  In 
Equation  74  can  be  obtained  by  multiplying  by  each  distribution  function 
and  integrating  over  (t)(y)  and  g . This  yields 


J OJ'P  <^(y)  - 'l>’(y)/2oJ  I d(j)Cy) 


^ Jexp  I -1  /l-R(C)*  go^  . g»/2  j dg  . 


(75) 


This  is  integrable  by  the  technique  shown  in  Appendix  A. 


R,3(5)  • exp  }.  oJ[l-H(e)]*/2  - 0»tl-R(5)*]/2  j 
■ exp  j-oJtl“R(£)l I . 


(76) 


(77) 


which  is  exactly  that  relationship  derived  by  Bromley 


(4) 


This  expression 


can  now  bo  used  In  Equation  71  to  obtain  the  angular  power  spectrum. 


3.4.2  Unmodified  Components 


The  expression  for  the  power  spectrum  cun  bo  written  in  a 
slightly  different  form 


P(0) 


lint 


x-«» 


A 

exp  j -o*tl"R(0]  - i 2iTeC/^  I d5  . (78) 
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The  approximation  will  be  used  that  the  autocorrelation  function  can  bo 
reasonably  represented  (reference  Equation  59)  by 


lUO  -•< 


1 If  e < Cl 


0 if  C > Cl 

where  Ci  • Insert  this  in  Equation  78. 


(79) 


pce) 


- ^^oxpl-  C'oJ,/2  -1  21T0C/X]  dC 

X 

* ':o»(2’t0C/X)  dC  (80) 


Cl 

>(0)  « exp  |^-C*oJ,/2  . i2nOC/x]  dC 

-Cl 


* Po  ii!  joxp(-aj)  lln-2n.e.x,/X.^:^sln  3it9,Cj/\  j 


lixpi'osslon  81  will  be  used  to  dotormino  the  unmodified  .sisnal  penetrating 
the  scattering  medium  and  the  scattered  signal  for  o,  » 1. 


3.4.3  Unmodified  Signal 

The  second  term  of  Expression  81  gives  the  unmodified  or  direct 
component  1’^.  Tor  0 ■ 0,  this  term  gives 


I 


P..  - 2Po  llni  (exp(-o*) 

( 9 A 


Xxw 


C82) 


m 


which  appraachos  tia  infinite  unguliii*  power  density.  Of  course,  the  Incident 
plane  wnve  had  an  infinite  anKular  power  density  there.  However,  If  wo 
perform  an  intosral  over  9 before  lettinn  x-+«  , we  will  obtain  the  power 
in  the  unmodified  beam.  (This  power  la  concentrated  in  an  angular  width 
of  the  ordcc  of  A/x  , which  accounts  for  the  constant  unscatterod  power 
flux  as  x-»«  .)  io  accomplish  this,  redefine  as 


llm 


fy  expi-oj)  JO 


The  direct  beam  will  be  defined  as  being  in  that  angular  width  Gj 
8ucl»  that  2ir0jX/X  >>  1 and  2irUjCi/^  1 ■ This  is  eqviivulent  to 
restricting  the  geometry  such  that  the  scattering  medium  is  at  least 
several  decorrelation  widths  wide.  With  this  restriction 


s^\  2ii0j^A  ,, 
■ |t0'  “ 


dO 


-0 : 


Pd  e 


18)) 


l«.8) 


The  last  term  In  liquation  84  is  negligible  based  iqion  tlio  restriction 
2tt0  Ci/i  1 . donseiiuontly,  the  direct  immodiflod  component  Is 
always  c.'  if  tiincs  tho  incident  power  and  can  he  roprosontod  liy 
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P„  - Po  0 S(0),  (86) 

where  6(0)  ia  u delta  function. 

3.4.4  Large  Phase  Perturbation  Scattered  Signal 

tf  the  phnao  perturbation  o^  ia  large 

and  in  Hquation  81,  the  liniita  of  integration  can  bo  extended  to  infinity 
with  little  error.  Integrating  thia  term  according  to  the  meana  of 
Appendix  A givoa 


Pi(0) 


O’ 


/^1t 


(88) 


a normal  d latrlhiition  with  variance  o^.  It  can  he  noted  that  aa  it  atanda, 
the  Integral  of  Pj  over  0 Is  Po,  the  total  Incident  power. 


The  term  involving  the  limit  on  x is  already  accounted  for  in 
the  iinscatterod  beam. 

The  final  term  in  liquation  81,  involving  ain  (SttGCj/^)  , can 
slightly  add  or  subtract  from  the  scattered  component  but  on  the  average 
subtracts  and,  indeed,  when  integrated  over  0 , removes  from  the  scattered 
beam  exactly  the  amount  of  power  that  is  In  the  unscattered  beam.  Its 
magnitude  la  always  leas  that  20  percent  of  the  main  scattered  component  and, 
for  valuta  of  In  excess  of  about  2,  la  loss  than  1 percent. 
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Consoi-iuontly,  little  error  Is  introduced,  colculations  are  simplified  and 
energy  is  conserved,  if  the  uvorogo  correction  due  to  the  lust  term  is  used 
for  nil  tingles.  Thut  is,  the  power  distribution  of  Eciuution  88  is  multi- 
plied by  1 - oxpC-o*)  . 

Using  this  lust  approximation,  the  total  power  angular  distribu- 
tion function  for  the  case  a,  » 1 can  be  written 


I <^^0  i 


3.4.5  Small  Phase  Variance  Scattered  Signal 


To  investigate  the  case  ot  small  phase  variance,  it  is  more 
instructive  to  return  to  the  basic  expression  for  tho  angular  power  spectrum 
Equation  78,  and  expand  tho  exponential! 


-oJU-lUOJ^  -oj 

3 » 0 + 


RIO  . 


When  tho  approximation  of  liquation  90  is  used  in  liquation  71, 
the  integral  of  the  firs*'  term  loads  to  tho  delta  function,  and 


! -0,  1 ‘1’  /*“’  1 
P(0)  - P ’o  6(0)  + ^ -i-T / R(C)  cos(2iiOC/A)  dC 

I ' ^ 


At  this  point  it  becomes  obvious  that  for  this  weak  scattering  case,  the 
angular  spectrum  is  not  necessarily  gausslun.  Thut  the  spectrum  was 


I 


guuasiun  for  strong  scattering  was  a result  of  o^[l'R(C)]  being  largo 
before  terms  in  11(5)  groutor  than  quadratic  were  Important.  I’or  weak 
scattering,  n broader  portion  of  the  autocorrelation  function  will  bo 
significant,  As  examples,  let  us  consider  a medium  of  uniform  gaussian 
rods  and  one  following  the  experimental  distribution  function  of  Equation 


i‘or  a medium  of  uniform  strlations  tho  autocorrelation  function 


11(5)  - exp[-5V4o*l  . 


(I’(cj)  in  Equation  42  is  6(0)  .)  This  can  be  expressed  in  terms  of  tho 
variance  of  phase  and  phase  gradient  by  making  use  of  the  relationship 
11"  (0)  » and  liquation  4?! 


II"  (Oj 


- . -IL  . 


IVlion  this  is  used  In  Expression  92  and  tlio  result  inserted  in  liquation 
91,  the  result  Is 


I ' ^'^“o  ( 


In  this  case  tho  scattering  does  happen  to  ho  gaussian,  but  with  an  angular 
variance  Increased  from  Op  to  • 

I'or  the  experimentally  derived  Chosnut  distribution,  the  auto- 
correlation function  is 


.If  «»«  «.»!■« 


This  cun  also  bo  Intogratod  explicitly  In  Gqiiutlon  91.  The  result  is 


P(0)  - P, 


(97) 


fn  this  case  tlui  distribution  Is  exponential,  not  giiusslan,  Init  there  Is 
little  power  in  the  scattered  component  because  oj-  is  assumed  to  ho  small 


I; 

P 
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For  weak  scattering,  the  energy  is  more  broadly  spread  than 
defined  by  the  parameter  Xo^,/2ir  , which  has  been  denoted  , The 
spread  is  increased  by  a factor  of  the  order  of  l/o^  , For  the  uniform 
gaussian  striatlon  model  the  variance  of  the  scattered  power  can  be  com- 
puted explicitly  for  al I . To  accomplish  this,  the  exponential  In 
Finuation  78  ciai  be  expanded  In  a Taylor  series  and  the  result  Integrated 
term  by  term.  Thus,  because  R((!:)  ■ exp(-?*/‘lo‘')  , 
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cosC2iieC/X)  d? 
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qS  = — i_  i . (104) 

Stt*  0*  1 - CXpC-CI^) 

Althoiinh  It  Is  clear  by  inrtiiectlon  of  laiimtion  1^9  that  even 
for  the  homoscnoous  striutlon  onvLronmont  the  power  JUtrihution  is  not 
i;au8slan  for  ».  I,  for  radar  modeling  purposes,  the  correct  char- 
acter to  the  signal  can  probahly  he  obtained  by  use  of  a gaussian 
distribution  with  a variance  given  by  lUiuation  104.  The  striatlon 
IKirainetcr  o can  be  taken  as  an  effective  size,  both  the  Impact  of 
the  distribution  on  the  riular  monopulsc  system  and  the  choice  of  an 
effective  size  will  be  treated  In  following  sections. 


3.4.7 


Summary 


One  way  of  choosing  an  effective  striation  size  Is  to  match  the 
autocorrelation  function  curvature  at  the  origin 

k'.'iO)  = - a - . (li 


if  this  is  Inserted  In  ihiuation  104,  the  result  Is 


Wrltlnj;  O*'-  i»\  this  I'oriti  cloutly  i I liisi  viitos  tliiit  t'or  tho  Boi'ui'Ul  o,aso  it 

has  tho  I'l'oiioi'  limiting  values  of  I'or  largo  a|  aiul  small 

0^  aiKl  t'or  tho  imiforiit  strlatlon  oaso  is  oovi’oct  throughout. 


It  Is  usoCul  to  think  ot'  tho  signal  as  ooiii|>rls lag  a true,  uiulls* 
tortoJ  signal  duo  to  tho  oniuodl fled  ooinpoiioat  and  a nolso  signal  duo  to  tho 
randomly  fliiotuat I ng  soattorod  oomponont.  In  those  terms,  tho  resulting 
slgnal-to-nolso  ratio  duo  to  tho  scattering  Is 


S/N  ^ 


oxin-i5|i 

oxpT-l^'f 


(108) 


3.5  Monopulse  Response  to  Scattering 

i'ho  important  i|Uost  Ion  to  answer  is  the  rosponso  ol'  a pioeo  of 
onulpment  to  the  signal  produced  hy  tho  scattering  onvi roniiioitt . It  will  bo 
assumed  that  tho  radar  pulse  Is  sufficiently  long  .that  tho  signals  propagated 
I'v  tho  various  scatter  paths  are  contomporarv  In  the  signal  processor.  The 
model  to  ho  developed  will  bo  th.it  of  a phase  comparison  monopulse  system. 

A "strong"  scattering  case  model  has  boon  previously  dovelopod, 
Tlie  analysis  was  specifically  applicable  under  conditions  satisfying  geo- 
metric optics,  ihe  analysis  hero  will  use  the  wave  front  Integral  and  thin 
pliase  screen  aiiproxlmat  Ion  that  Is  appllcal>le  under  the  single,  loss  strin- 
gent condition  that  tho  phase  shift  variance  is  largo.  (The  results  will  bo 
seen  to  be  Identical  to  those  of  the  earlier,  more  restrictive  derivation.) 
The  analj'sis  will  then  bo  extended  to  conditions  of  small  phase  shift 
variance, 


39 


! 


3- 5.1  Sum  and  Dlffe^rcnce  Channel  Signal  Formulations 


IMlvuro  2 hIvowh  the  no,lM#cti'y  ot'.tho  cunFiKiiViit  ion.  The 
nntonna  Is  h plmscd  uvniy  point Vnj!  on  nn^jlo  a fpom  thd’ sourco.  'I'hc 
total  phase  shift  roKitlvo  to  the  direct  path  Ls 


^ ^ sin  a j + 


|10«V 


T \ r 


Srt  j 


+ (}>, 


(Ill)) 


wluM'c  (1  tonii  In  t*/Ar  has  hcon  dropped  hceimse  It  Is  small  In  the  fur 
tone  of  the  iintoin\iii  sin  o.  has  hcon  reprosoiuod  t>y  u (a  wrltlnp 


simplification,  not  really  a raathomat  lea  1 approximation);  h"'  “ iC’  * >'’* 


'1' 


as  hoforci  and  hlRhor  order  terms  In  .;/r  , y/r  , and  y/H.j,  have  been 
dropped  lu'Ciiuse  scattcrlnj;  anj'les  have  been  assumed  to  Ih;  small. 


Figure  2.  Antenna  - scattering  configuration. 
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l^riinir  n'  'iilSir  i . .1  f Jl  . 


Tor  an  nntonna  of  width  d , tht*  nomilizod  iiionopulao  humi  channel 
ainplttado  1? 


z-r 

r«  111).,. 

dy  ds 

(111) 

•'z— d/2 
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•'-d/2 

/"oxpl- 
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-a)  z + UM  dy  dz 

(112) 
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(113) 

1 "d  / 

r"'  sin  Y f • 

L 

j — • expl 

XU 

1(1)  1 dy 

(IM) 

To  make  the  iiiiitlu'iiuit les  traetahU'  tlie  antonnii  pattern  l,»ln  x/x)  will  be 
approx Imated  by  a paasslan.  In  the  I’lital  application,  some  of  the  jjausslan 
factors  cancel,  so  that  the  approximation  Is  only  einilvalent  to  flttlnn  the 
antenna  pattern  with  a gausslan  locally  about  the  ri'ulon  of  maximum  return. 
Us  lap 


I 


i 


UlRl  a 

ij 
5 

liipiatlon  IM  becomes  i 


In  actual  practice,  the  pausslan  may  be  a better  approximation  to 
an  antenna  pattern  because  feed  taper inp  is  used  to  suppress  the  side  lobes 
associated  with  uniform  1 1 lumlna.t ion  of  a rectanpular  aperture, 
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I'or  future  writing  simplicity  denote 


b - ft/XK 
c * Tr*d*/6X*. 

The  difference  channel  output  will  ho 
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- y*  ( ^ - a)  expl-  iby*  • | C ^ - a)*  + dy 
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whore  u gausslan  approximation  to  the  antenna  pattern  has  been  used  once 


3.5.2  Sum  Channel  Power 

For  normalization  purpoaea,  compute  the  undlaturbed  aum  channel  vol- 
tage and  voltage  variance,  which  ia  proportional  to  the  puwer,  for  0 ■ 0, 
a ■ 0 , 


■ d J oxpt-lby*  - cy*/r*  ] dy 


Po  ■ liSol* 


E.i;.. 


36X‘'r''  X'R* 


In  tho  nntonnu  fur  field  d^/Xr  « 1;  hence  Equation  124  is  nearly 

Po  ■ d*AR  • (1 

Consider  next  the  mean  value  of  the  sum  channel  voltage.  The 
distribution  function  of  tho  phase  shift  of  tho  phase  screen  is 


so  thut 


/Iv  0 


^ — oxp(-  (fi*/2o’) 


" f f exp  j - iby*  - c(  ^ - o)*  + 1<|>  - i|)V2o?|  dy  d()) 

J J a ^ ^ I*  I 


• 00  aOO  |h 


This  Integral  is  the  form  of  that  of  Appendix  A.  Because  tho  «'orinblea 
are  independent i integrate  on  first.  Thl.s  loads  to 

n*a*d* 

f . / -t  / 

A 

wl\cre  the  far  :unc  approximation  (c  « Ihr)  has  l^ocn  used. 

Caro  must  he  uxorclsod  either  In  tho  Interpretation  of  t or 
because  of  antenna  filtering  effects.  If  a*  is  based  upon  nil  spatial 
wavelengths  that  might  ho  present  in  the  environmental  structure  Including 
those  compnriihlo  to  tlic  subtense  of  tho  antenna  beam  (i.o.,  (1  varies 
systematically  over  the  antenna  heam) , there  Is  a portion  of  tho  signal  pltase 
variability  that  Is  constant  over  the  beam  at  a particular  time  but  varies 
with  time.  A phasor  plot  of  tho  signal  may  just  bo  a nearly  constant  ampli- 
tude vector  rotatittg,  possibly  slowly,  about  tho  origin.  In  this  case,  It 
is  tho  vector  average  of  the  signal  which  has  the  value  given  by  i.c|uatlon  12H. 
liven  a nonscintJ  Hating  signal  may  have  a nearly  aero  moan  due  to  a slowly 
varying  gross  phase  shift. 

A possibly  more  useful  approach  is  to  Include  a spatial  wavelength 
cutoff  In  the  computation  of  such  that  the  ensemble  average  of  lU|uatlon 
127  has  some  valid  moaning  in  the  context  of  a single  measurement.  Tho 
longer  wavelength  structure  would  then  he  accounted  for  by  means  of  a 
moan  phase  shift,  gross  refractive  angular  error  and  focusing  gain  (to 
roprosont  mean  phase  shift,  moan  first  spatial  derivative  and  mean  second 
derivative).  Tho  proper  way  to  separate  the  effects  that  arc  duo  to  random 
variations  within  u boamwidth  and  those  due  to  correlated  variations  hus 
not  boon  developed, 
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Similar,  but  even  more  stringent  limitations  may  apply  in  the  case 
of  small  phase  perturbations  or  small  mean  scattering  angle  because  the 
region  of  significance  may  be  even  smaller  than  the  antenna  beamwidth. 

It  may  be  controlled  by  the  mean  scattering  angle  or  a Fresnel  tone. 

If  the  one-way  power  gain  of  the  antenna  relative  to  boresight  is 
called  liquation  128  can  also  be  written 


— -4 

d /WT  e ^ ClCa)''*  . 
The  square  of  the  real  part  of  S is 


C(a3 

which  (by  comparison  with  liquation  85)  Is  the  dli'oct  uniiKHllfled  power 
attenuated  by  the  off-axis  rejection  of  the  antenna, 


Compute  next  the  mean  power  received 


1:1: 


d*  y*  J'  exp  jib(-y®+z®)  - c(  ^ - a)*  - e.(  - a)® 


+ iijiCy)  - i(j>(a)  dy  dz  (1.51) 


Let  a » y + 5 


^ ♦ !(<!>(>’)- Ky  ♦<:))}  Uy  d5  . (i: 

The  only  fiictor  thut  nocUa  uvoriiRlnfi  over  the  ensemhlo  is  that  Involving 
the  acx'oen  phuao  shift  a.  This  was  previously  dofinod  CHi|uution  73)  as 
thu  autocorrelation  function  of  the  electric  field  K|,Ct)  and  evaluated 
to  bo 


« cxp{-oJ[l  -IKOlj. 
Uquation  13 J can  then  bo  rewritten  us 


ci: 


noo 

cxpk'lib--^]  ♦ i-ll^by-.^(^-cx)] 

W ' 

- «)=*  -ojl  l - lUfJj}  dy  dC  . (I, 


In  the  limit  thut  » 1 , tho  integrand  is  sharply  petikod 
about  C H 0 where  tho  term  l*<  nearly  zero,  llonseiiuently , 

the  Taylor  scries  approximation  can  bo  used  for  lUC)  (see  liquation  59) 
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With  this  npproxiinntion,  Kquntlon  134  enn  he  Integrated  first 
over  C then  >•  to  produce 


PjL 


y.;^”'Ts(er..iiKi 
^ 0 


It*  a* 


(136) 


where  terms  small  In  the  far  zone  hnvc  been  dropped,  Og  ■ 
and  00  « X/d  (the  position  of  the  first  null  in  a rectongular  antenna's 
side  lobe  pattern).  This  can  also  be  written  in  terms  of  the  antenna  gain 
at  angle  ct 
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(137) 


wltero 


u ■ 4tt*  0^/6  0|, 


(138) 


It  can  bo  shown  that  if  Ti*a*/30?  Is  greater  than  0,5,  scattering 
limy  increase  the  signal  strength  over  the  value  ohtoinod  in  a non -structured 
environment.  Unsically,  this  moans  that  if  the  target  is  sufficiently  far 
from  the  beam  pointing  direction,  scattering  may  provide  a higher  gain  patl\ 
for  the  energy  than  the  direct  path. 


It  is  clear  that  the  mean  power  greatly  exceeds  the  power  in  the 
unmodified  component  if  o^»  1 . Furthermore,  Inspection  of  Equation  116 
shows  that  the  signol  is  effectively  composed  of  many  randomly  phased  com- 
ponents because,  for  the  assumptions  used,  the  li)>  term  can  vary  over  many 
radians  within  the  boamwldth  and  Fresnel  zones.  As  a consequence  of  the 
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ccritrul  limit  theorem,  In  thlii  limit,  the  »uin  channel  voltauo  will  have  In- 
phuiio  unJ  (.(uudraturo  cor\poi\cnta  that  are  ouch  Independently  ){aua!tlnnly  dis- 
tributed with  p variunco  equivalent  to  half  the  received  power. 

rrr  (j*  « 1 , ju  alternate  assumption  allows  intciiratlon  of 
Equation  l.M.  The  exponential  factor  Involvlnn  the  autocorrelation  function 
cun  be  expanded  us  u Taylor  sorios  In  with  the  retention  of  only  the 
linear  term  In  o?  . 1‘hus, 

-oj  »j«i« 

0 ^ ■ e ^ e (1311) 

- e |l  lUCj)  . 040) 

The  autocorroliit  Ion  function  can  then  be  approximated  l>y  a yaussian, 
oxpi-af,*/2)  , and  tl\e  result  Inserted  In  liquation  ISd.  The  naussliin  auto- 
correlation function  is  exact  only  for  uniform  site  strlatlons  if  the 
profile  Is  kiaussian  (this  can  be  deduced  from  liquation  -Id  and  will  be 
discussed  later).  However,  Inasmuch  us  It  Is  the  only  pliyslcally  roason- 


able  function  that  renders  liquation  ISd  Inteurablo,  wo  will  use  it  in  the 
I hoi)es  of  obtuininji  at  least  a reasonable  estimate  of  the  scattering, 


This  Is  of  the  same  form  as  previous  integrals.  If  we  make  the  fur  zone 
approximation  (c*  « b^r**),  match  the  curvature  of  the  autocorrelation 
function  of  the  phase  at  the  origin  (a  » tl'®  limit  of 
small  , approximate  o^  by  oxpCu^)-li  wo  obtain  for  the  integral 


exj:) 


a* 


L .U*|l  + 


1 


b r*  0* 


(142) 

This  can  also  be  written  in  terms  of  the  antenna  beamwidth  0^  and  mean 
apparent  scotterina  variance  a*. 


(143) 

By  roforence  to  liquation  130,  it  cun  bo  seen  that  the  first  term 
is  exactly  tho  power  associated  with  the  unmodified  component,  and  the 
second  term  must  be  that  produced  by  the  random  noisy  component.  If  tho 
derived  value  of  tho  moan  scattering  variance  is  obtained  from  liquation  107 


1 - oxp  (-0*) 


(144) 

(145) 
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und  tho  asymptotic  forms  for  large  and  amall  are  used  in  the  expressions 
for  strong  scattering,  Equation  136,  and  weak  scattering,  Hquation  143,  the 
result  in  both  cases  is 


P « Pfl 


4. 

/l 

V 3 G? 


exp 


(.146) 

Thus, we  have  a single  expression  valid  for  both  strong  and  weak  scattering. 


3.5.3  Sum  Channel  Amplitude  Distribution 


Por  nti'ong  scattering,  it  was  shown  that  the  sum  channel  signal 
was  composed  of  two  independent  gausslnn  variates  In  quadrature.  Mow  Is 
It  fur  weak  scattering?  P.xpand  tho  basic  sum  channel  voltage  given  in 
P.quntlon  116  for  small  phase  shifts. 


on 

d J" oxpj-lby*  - c(^  - a|*|  |l+li>  ♦ ...j  dy  (147) 


w 

1^0  d ycixpj“iby*  - c(^  - aj*  j H dy  (148) 


Tl\is  shows  the  sum  signal  is  nearly  tho  undistortod  signal  plus  a weighted 
integral  of  tho  phase  shift.  The  integral  is  tho  limit  of  a sum.  Conse- 
quently, tho  integral  ovor  tho  weighted  gaussian  variate  will  be  n 
guussian  vnrlato.  Consider  the  integral  as  u summation.  I'hore  will  bo 
regions  of  y in  which  tho  phase  is  relatively  constant.  Call  the  width 
of  such  B region  Cj  and  take  it  to  bo  the  distance  for  tho  autocorrela- 
tion function  to  drop  to  0.5, 
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The  character  of  the  integral  of  Equation  148  ij  that  there  is 
general  attenuation  near  y-0  of  exp(-ca*3,  which  is  the  antonna  response 
in  the  target-source  direction, and  a continuously  changing  phase.  The 
regions  of  coherent  plasma  phase  shift  will  be  coherent  in  the  Integration 
only  if  the  spatial  phase  change  across  the  region  5i  is  small.  Outside 
that  region,  the  phase  becomes  randomized  and  does  not  add  coherently. 

Tile  phase  shift  gradient  due  to  )ihysleal  path  length  Is 


for  estimation  purposes,  require  that 


d(^^ 

W 


Ki 


(1513 


This  establishes  n maximum  region  of  the  phase  screen  y^^  outside  of  which 
the  spatial  variation  of  phase  is  so  rapid  as  to  prevent  coherent  phase 
addition.  Thus, 


^m  “ 


(152; 


If  the  structure  dimensions  are  small  compared  to  the  first 
Fresnel  zone,  the  total  spatially  induced  phase  shifts  involved  in  y^^^ 
are  largo  compared  to  2ti.  The  contribution  of  the  phase  shift  should  he 
iioaii;.’  random  so  long  as  the  total  pluiso  sliift  in  y^^^  is  large  compared 
to  2tt,  Across  each  element  ^i,  the  contribution  to  the  integral  of 
liquation  14S  has  a variance 


i 

■4 


AP  ' d^  Cf  i;(a) 


(1.5.-; 


of  this,  because  the  phase  shift  associated  with  the  exponential  factor  in 
liiiuution  148  spans  a lurjio  variation,  on  the  average,  half  the  contribution 
will  bo  in  phase  with  the  unmodified  component  and  half  in  quadrature. 


The  number  of  sucli  Independent  sampler,  will  he 


11  o 


/f,, 

' m • 1 


(154) 


The  variance  of  the  sum  of  n samples  is  n times  the  variance  of  the 
single  sample.  Thus,  the  In-phuse  component  of  scattered  power  is 


1',^  ™ I d' Mlojllu*)  (15S) 

» |l'ooJ(Hot)  . (,15b) 

The  result  of  liquation  iSo  is  comparable  to  the  scattered  component  in 
liquation  14(1.  It  differs  in  that  the  sum  of  the  in-phase  and  an  equal 
magnitude  quadrature  power  Is  half  of  that  obtained  in  the  detail  calcula- 
tion, which  Is  reasonably  close  considering  the  hoc  assumptions  used  In 
the  ertide  derivation,  It  also  differs  by  the  neglect  of  —s—\lr  compared  to 
Oj  , but  this  is  valid  for  • If  plasma  structural  dimensions  are 

small  compared  to  the  first  I'rosnel  zone. 


The  inforiimtion  added  by  this  quantitatively  more  crude  analysis 
is  that  the  scattered  component  consists  of  independent  in-phase  and 
quadrature  gausslun  components.  This  of  course  leads  to  Kleiun  statistics 
f""  the  power. 

In  addition  to  the  evenly  divided  components  tliorc  Is  an  addi- 
tional contribution  of  a random  component  in  quadrature  to  the  direct  signal 
tliat  is  duo  to  the  largo  extent  of  the  first  ircsnel  tone,  liowever,  it  cun  be 
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shown  to  be  smrillor  by  the  rntlo  of  the  plasiim  .structural  Jimension  to  the 
first  bresnol  zone  size.  This  has  been  reforrcd  to  ns  I'resnel  filtering 
and  may  bo  significant  for  some  coiablnutions  of  and  spatial  autocorre- 
lation function. 

3 • 5 • 4 Difference  Channel  Power 

Consider  now  the  difference  channel,  for  which  the  signal  voltage 
is  given  by  lUiuatlon  l«’l.  If  the  phase  perturbation  Is  zero,  the  undis- 
turbed reference  signal  is  obtained.  Thus, 

Ao  £ (^-  j o.xp  y'(-lb  - + .^'-y  - ^..a’  dy  ll.‘l7i 

\ 

la  the  far  zone  of  the  antenna,  this  Is 


It  can  be  seen  by  reference  to  lUpiat  Ion  1J.1  that  the  difference  channel  is 
in  quadrature  to  the  sum  channel,  l■urthern^ore,  tlie  exponential  term  Is 
just  the  antenna  gain  of  each  half  the  antenna  used  to  form  the  difference 
hoani,  or  si)eelf Lealiy  fITCa/Jj  . The  other  term  In  the  parenthesis  is 
Just  the  small  angle  approximation  to  the  noriiial  difference  channel 
function  sln(pu)  where  p is  iid/2.\  . Thus, 


AjA*  = I’,  sin^tpa)  iHa/Jj  . (.151 

In  the  presence  of  phase  fluctuations,  the  variance  of  the  dif- 
ference channel  voltage  will  differ  from  that  given  hy  liquation  1.11),  It 
can  he  comimtod  In  a maimer  similar  to  that  of  the  sum  channel: 


....... 


■ f [(r  *'^f  (r  * “)*]*  I '^y  ‘'’■ 

(IftO) 

If  titc  rtiihrttltutlon  : " y+C  nscil  iiiul  the  phiiso  iliffercnco  uvoriigcii 
over  the  oiiscmblc,  the  result  Is 

’ X^//  I)  } (i”  - A^r)  ■ fJ- 


- I n*  - U-iUOl  Jy  tl?: 


whore  ii  • ^.-  - a. 


I'or  the  eiise  of  liu’ne  phoso  pert urbiit ions , the  term  involving 
the  autocorreliit  Ion  function  can  he  npproxiimiteil  hy 


«JU  -lUOl 


r?  2 0'^ 

f.  > . l/> 


As  was  dlscusscil  for  the  sum  channel  power,  the  Intenral  on  C in  lii|uatlon 

U)1  can  lie  broken  at  f,*  « Ixtenral  for  laree  just  nlves 

the  uninotlLfied  component  times  nn  attenuation  factor  0''^$  plus  an  attenu- 

-ci* 

at  Ion  of  the  scattereil  component  that  avernyes  l«e  'll, 

bet  us  concentrate  on  the  portion  of  the  Integral  near  tlie  origin. 
As  a mutter  of  fact,  the  term  exp(.-o* liecomes  so  siiiall  outaUle  the 
Interval  of  concern  that  the  limits  on  the  integral  can  he  extoncloil  to 
Infinity  with  little  Introiluctlon  of  error.  Thus,  the  mean  scattered 
component  is 


"’y  jf  I )®’‘l’  { [4^  " -f^  • ^1^] 


+ C 


[i.^by  ‘ S]  * J>’ 
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Straightforward  application  of  Appendix  A to  first  C then  y yields 

'’"!(i5&)  * t(  ?)! 


23^ 

s 


('  ‘ r 


'I  1 

0x11  j 


TT*  a* 


uei  (1.  J' 


Clf.4) 


The  first  term  In  the  curly  brackets  replaces  what  was  the  diffcronco  signal 
sli\*(pa)  and  comprises  a component  that  Is  a function  of  the  direction  off 
the  inonopulse  null  and  a component  dcpeiulont  only  on  the  scattering,  which 
makes  the  object  have  tl\e  appearance  of  moving  around  off-boreslght. 


it  is  interostlng  to  rewrite  liquation  164  in  terms  of  an  effective 
angle  «'  ■ a / /l'+u7'f • As  will  be  seen  the  effect  of  the  scattering  can  bo 
interpreted  as  moving  the  apparent  position  of  the  source  closer  to  borcsight 
by  the  factor  l/ZT+rrA.  liquation  164  becomes 


A'A'*' 

s 


i,j  c;(a'/2,)  slnAl'rt'. 


I*  5Tr 


U65) 


This  can  lu’  interpreted  In  terms  of  a power  loss  due  to  the  scattering  and 
the  reduction  in  effective  boroslght  angle,  and  shows  that  the  relative 
magnitude  of  the  iion-horoslght  dcpeiulont  term  is  Just  the  ratio  of  the 
variance  of  the  phase  screen  scattering,  as  viewed  from  the  radar,  to  the 
square  of  the  apparent  off- horeslght  angle. 
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In  8umumry,  for  o*  » 1,  the  difference  chunnel  sliinnl  comprisofi 
an  attenuated  unmodified  component  Ao  oxp(-o^/2)  and  ti  scattered  random 
component  that  can  he  argued,  in  n manner  siinilnr  to  that  for  the  sum 
channel,  to  have  Independent  random  in-phase  and  quadrature  components 
with  variance  (1  -oxpC-o^)]  A'A"*/2  . 

Consider  now  the  weak  perturbatii'ii  case,  As  for  the  case  of  the 
sum  channel,  the  random  pha.so  effects  will  he  approximated  hy  taking 

o.xp  j -o|[l-IUfJ]  I " oxp(-o|5  + [l-oxp(-o*i]  fU(C)] , C16(0 

With  tlio  use  of  this  approximation,  liquation  Ihl  hecomes  the  sum  of  two 
Integrals.  The  first  is  exactly  expi-ap  AoAo> 

If  the  autocorrelation  function  Is  once  more  approximated  hy  a 
gaussiaii,  U(C)  • exp(,-o*,C*/2o|) , one  sees  that  the  second  Integral  Is 
Identical  in  form  to  that  integrated  for  strong  scattering  except  that 
wherever  appeared  originally,  aiipears,  As  was  dis- 

cussed earlier,  this  is  equivalent  to  replacing  hy  ip  wliere 


Wltli  this  modification,  the  same  equations  apply  in  both  limits  a*  c<  i 

and  0?  » 1. 

'll 


3.5.5  Sunwiary  of  Monopulse  System  Output 

In  summary,  this  gives  a unified  moans  of  computing  the  sum  and 
difference  channel  signals.  The  equations  are  explicitly  valid  for  largo 
and  smail  values  of  o.  and  It  Is  assumed  that  thev  can't  he  too  far  off 
for  a.  % 1.  The  environment  is  defined  In  terms  of  and  either  o*;, 

ifi  (|)  9 ' 


or  0^  which  are  related  by  ■ IUo^,/2iTr.  Means  for  computing  these 
parumotors  from  un  environmental  doscriptlon  were  given  in  an  enrlier 
section  and  specific  examples  follow  in  the  next  section. 

The  radar  system  Is  character  lied  by  antenna  boamwidths 
6j  ("l/d)  , the  first  null  in  a rectangular  aperture's  pattern,  and  0^^  , the 
• 3 dB  power  henmwidth  which,  for  a rectangular  aperture.  Is  related  to  0^ 
by  ■ 2.8Qo/t.  The  sum  beam  relative  power  gain  an  angle  a from  beam 
center  is  C!(0().  The  on-axis,  undisturbed  sum  beam  signal  power  is  Po  • 

The  undisturbed  difference  beam  power  is  Pc  n(ct/2)  sln*(pci)  whore 
p ■ TTd/2X. 


A parameter  characterising  the  relative  importance  of  scattering 
and  untoniui  rejection  is 


u « 


i!:JI 

30  0 


8.23  iF 


(UiH) 


Using  these  basic  Inputs,  the  sum  channel  .signal  comprises  an 
unmodified  component  of  power 


P^,  " Po  e.\i»  1-0*1  (Hal 


(109) 


and  two  orthogonal,  random,  independent  components,  onch  guussianly 
distributed  of  variance 


(l-oxp(' 


o*)l  oxp 


2 . Qg* 


u 

■Hu 


/O'*! 


(170) 


The  difference  channel  comprises  an  unmodified  component  in  ciuadrature 
to  the  unmodified  sum  channel  signal  and  a scattered  signal.  The  un- 
modified component  has  a power  associated  with  It  of 
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AoAo  ■ I’o  oJ<p(-a?)  ri(a/23  9in*(px). 


(171) 


The  scattered  component  comprises  phase  orthogonal,  indopondont  random 
gniissitm  components  with  variance 


•>' " TjTTHp*  ^ * T t) i • 


It  should  ho  pointed  out  that  on  the  current  effort  possible 
correlation  between  the  sum  and  difference  channel  signals  has  not  been 
analyzed  and  the  difference  channel  analysis  has  assumed  that  the  target 
is  displaced  along  one  of  the  axes  of  the  monopulse  system  and  that  tho 
strlation  axes  arc  aligned  with  tho  other.  Analysis  to  overcome  these 
shortcomings  is  straightforward  but  there  was  insufficient  time  to  fully 
explore  all  aspects  of  the  processor. 

4.  SCATTERING  BY  VARIOUS  PLASMA  STRUCTURES 

This  section  will  Investigate  tho  relative  scattering  pi’oduced 
by  different  size  distributions  of  strlatlons  to  determine  the  sensitivity 
of  scattering  calculations  to  details  of  tho  distribution.  Tho  problem 
is  tl\ut  tho  distribution  Is  not  well  established,  so  if  there  were  a great 
sonsitivity  to  plasma  parameters,  there  could  be  Tittle  rclinnco  in  tho 
quantification  of  the  scattering. 


4.1  Model  Descriptions 

Four  models  of  tho  distribution  of  strlation  sizes  will  bo  con- 
sidered us  illustrative  distributions.  Tbc  particular  models  were  chosen 
for  simplicity  or  prior  uso  in  tho  literaturo.  Tho  models  are  single, 
uniform  size,  a modification  of  Chosnut's  distribution,  and  power  law  dis- 
tributions with  either  fixed  or  sizo-depondont  axial  electron  content. 
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4.1.1  Model  1 

In  this  distribution  all  strintions  have  the  same  characteristic 
radius  and  peak  oloctron  concentration.  The  probability  distribution 
function  of  the  I’ndll  is  a delta  function,  This  is  an  idenllied 

model  that  is  easy  to  visualize  and  to  work  with.  Because  it  is  so  simple, 
the  results  of  analysis  of  the  scntterinR  properties  of  other  distributions 
will  ho  expressed  in  terms  of  an  effective  striation  size.  That  is,  the 
characteristic  dimension  of  uniform  strlutlons  needed  to  produce  the  same 
scatterlnji  will  he  determined. 

4.1.2  Model  2 

This  distribution  was  inferred  by  Walter  Chesnut  from  nn  analysis 
of  photographic  data.  The  probability  density  function  is 

'■w)  aw 

Rcasonalilc  distributions  (i.o.,  finite  phase  shift)  imply  that  n exceeds 
4.  The  normalization  factor  a is 

pi-: 

M 

a “ ” — 

/ft' 


The  specific  power  favored  hy  Chesnut  was  n»li  so  that  a was  8/S/n. 
A plot  of  this  Function  for  n*.'!,  b,  7 and  8 is  shown  in 

I'Ltiiirc  .1. 


n-l 


^)l 


[n-S) 


, n even 


F ■ 


n odd. 


(174) 
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To  relate  this  fuiutlon  to  the  plusiiiu  Instability  analysis,  which 
yields  u minimum  cutofT  wavelength,  it  is  necessary  to  develoi>  a definition 
of  the  ohurueteristlc  striution  size  . if  the  cutoff  size  is  taken  to 
ho  that  hiivinn  a tenth  the  probability  density  at  the  peak  of  the  spectrum, 

Oj  - 1.8C,  o^n''"  (170) 

Js  a reasonable  fit  for  the  data  of  I'iKurc  .1, 

As  was  pointed  out  earlier,  this  spectrum  of  radii  applies 
specifically  to  fireball  roRlons,  some  of  the  observed  spatial  power 
spectra  were  best  matched  with  the  sum  of  two  distributions  of  the  form  of 
liquation  173,  and  the  analysis  was  not  very  precise  about  the  exact  value 
of  the  exponent  of  a^/o  — a value  of  (>  seemed  to  fit  the  data  us  well  as 
any.  However,  it  is  a convenient  Intexrable  form  tliat  fits  the  data 
reasonably  well,  Also,  because  the  phase  shift  approaches  a yausslan 
random  sequence,  the  variances  of  phase  shift  and  of  derivatives  of  phase 
shift  for  a distribution  that  Is  the  sum  of  several  components  arc  Just 
the  sum  of  the  variances  for  each  of  the  several  components. 

Another  way  of  viewing  this  distribution  Is  ihut  It  is  a power 
law  distribution  with  the  short  spatial  wavelength  end  cut  off  by  the 
exponential  rather  than  by  an  abrupt  limit.  As  such,  It  is  probably  more 
realistic  than  the  pure  power  law  distribution,  i'lio  continuous  function 
has  the  advantage  of  allowing  the  limits  on  integrals  determining  plasma 
scattering  parameters  to  be  o ■ 0 or  ® , 

The  axial  electron  concentration  Is  assumed  to  be  constant, 
independent  of  striution  size. 


4.1.3  Model  3 


This  1r  a power  law  Uiittrlbutlon  with  n normalized  distribution 
function  of 


I'lo) 


1 


(1771 


(s  - 1)  „ IS  - 1) 


for  cj,  < a < and  zero  outside  those  limits.  This  reflects  tho  situa- 
tion that  below  a minimum  cutoff  size,  plasma  perturbations  are  stable  and 
will  not  urow  ami  that  structures  larger  than  radar  beam  dimensions  or 
plasma  gradient  distances  have  either  no  Influence  on  the  scattering  prob- 
lem or  no  moaning  in  terms  of  the  plasma  description  used  here. 

There  Is  little  reason  for  u particular  choice  of  power,  lixcopt 
for  the  arbitrary  cutoff  for  large  structures,  this  distribution  would 
also  produce  infinite  phase  shift  variance  for  s£4. 

Tho  axial  electron  content  is  taken  to  bo  uniform. 


f 4.1.4  Model  4 


I 

I 


Tills  model  has  the  same  distribution  function  us  Model  3,  but 
tho  axial  electron  concentration  Is  taken  proportional  to  tho  strlation 
size. 


n,  ■ nj.0  . C178) 

This  is  based  upon  an  intuitive  feeling  that  any  growing  perturbations 
should  pinch  off  and  stop  growing  once  the  perturbation  amplitude  becomes 
a ronsonuble  fraction  of  a structural  wavelength. 


l 
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This  was  the  model  used  In  the  initial  IIARC  analysis.  If  all 
possible  perturbation  models  are  allowed  in  a cylindrlcally  posed  problem 
and  modes  are  allowed  to  decorrolato  in  a wavelength,  a value  ot*  four  is 
deduced  for  s . 

4.2  Plasma  Parameter's 

It  will  be  assumed  that  in  reasonably  homogeneous  regions  the 
local  electron  concentration  variance  o*  and  u cutoff  wavelength  are 
available  from  either  plasma  analysis  or  experimental  data.  Those  will 
bo  used  to  compute  the  variance  of  phase  shift  and  phase  shift  transverse 
gradient  which  cun  then  ho  integrated  along  sight  paths, 

The  result  cun  be  expressed  in  terms  of  nn  effective  uniform 
distribution  strlatlon  site.  The  pertinent  equations  are  repented  here 
from  previous  sections; 
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and  a combination  of  liquation  47  and  48 
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Model  1 


If  liquations  18.3  and  IHS  are  ova  luatod , It  Is  appai'ont  the 
strlatlon  size  of  the  uniform  strlatlons  Is  iUontloally  the  same  as  c 
and 

si  so  0 


0|.,  This  is,  of  course,  the  reason  for  having  labelled  the  uniform 


Two  Interesting  points  are  doveloiu'd  when  the  scattering  para- 
meters uro  exproasod  this  way.  The  concentration  of  atriations  and  the 
axial  electron  content  no  longer  appear  explicitly  in  liquation  182  or  184 
and  the  ratio  of  the  variance  of  phase  and  jihase  pradlent  Is  simply  related 
to  the  strlatlon  dimension,  liquation  182  also  Implies  that  the  variance 
of  the  scattering  angle  Is 
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4.2.2  Model  2 


The  effective  .striution  size  Is 


ou 

f •’  1 1*- 1 

0 

00 

f x”'"*  cxp(-x^3  ilx 



j x"*'^  exp(-x^1  dx 


where  x « Oo/o.  I'or  n even, 


c c ,n-o 


or,  for  n odd. 


These  lire  shown  in  I'ijiiue  T for  values  of  n from  4 to  i).  Similarly, 


oi  = I . «(i^  o''  r 


J x"'’’’  cxp(-x^)  dx 


J'  x""''  expt-x^}  dx 


Valni's  of  0 1.  are  also  shown  in  i'ljiiire  4. 


7 8 9 10 

Power  of  distribution  function  - n or  s 

Figure  4.  Effective  strlatlon  sizes  for  the  various  scattering  models, 
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4.2.3  Model  3 


TIk'  l.iito};i’iils  I'oi'  Moilol  3 nrc  slmpLc  iiiul  yiolil 


0^,  = — ,,,  . 

s-3  1-n  ■ 

S-'l 
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wIkM'O  11  Oi/O;;  — tlK’  I’lll'ln  I'T  tiu'  loKl'l’  illlll  lllHIOV  Ullt-ofl’  StrllltLo'.l  !«ir.l'Sl, 
1^11'  illstrllHition;i  wllli  n oiiiial  to  »,  3 or  •!,  the  rollowln}'  limit  ran  'm 
ui'Ovl  111  tlioso  osprossions 
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4.2.4  Model  4 

lor  Model  '1  the  otTort  1 vo  siioa  are 
11- 

0 r,  0 1 , lli'l'l 

3-s  l-n 


, , s-a 

(I'S  l-l) 

Tlioso  are  also  plot  toil  In  I'iiniro  •!. 

4.3  Comparison  of  Model  Results 

I’ornsal  of  the  data  In  l■l(;llro  I shows  somo  not  too  surprising 
troiuls.  Tho  larnor  strlatlons  are  ompliasrrod  (.with  rosultlnp,  larnor 
orfootlvo  slro.)  by  lower  Inverse  powers,  by  liaviii);  larger  axial  eleetroii 
eontoiU  iMoilol  'll  aiul  by  the  model  that  rolls  off  more  slowly  at  the  short 


wavo  ciitoft".  llowovcr,  it  is  also  noticeable  that  ti  lot  of  the  effective 
strintion  sizes  lie  between  1.5  and  4 times  the  cutoff  size,  Consoquently, 
if  un  effective  striatlon  size  i.s  taken  to  be 

- 0^.  ■ 2.. 5 (198) 

rousoiuiblo  results  arc  anticipalod.  This  value  Is  recommended  because  It 
is  near  that  derived  for  the  one  exiierlmental  distribution,  Also,  tlie 
iiumnotlc  field  should  continue  to  break  up  the  iurtier  striatlons  so  that 
low  powers  In  the  distrllHit ion  function  shouldn't  persist,  I'u. thermore , 
the  sharp  cutoff  that  leads  to  small  values  of  seems  physically  un- 
realistic, It  would  .seem  surprising  If  the  proj'er  effective  size  differed 
from  2,.S  0 by  a factor  of  2 and  It  most  likely  lies  within  a factor  of 
I.,";. 

Certainly  with  our  current  state  of  knowledge,  there  Is  no 
reason  for  iisslgnlng  different  values  to  and  Oj.  . 

It  should  he  remembered  that  o';,  and  a?,  vary  Inversely  with 
, Consequently,  the  scattering  uncertninty  introduced  by  the  striatlon 
uncertainty  Is  rodviced  by  a square  root  process, 

'fhe  pliase  variance  oj  v;iries  as  the  product  2o^,  oj , , Con- 
sequently, If  the  effective  striatlon  size  were  actuiilly  smaller  than  used 
here,  would  be  greater  and  would  be  smaller,  l)oth  l)y  aliout  the 

fractional  error  in  o . 


5. 


NUMMARY  AND  UTILIZATION  OF  FORMULATIONS 


It  will  bo  usaumorl  horo  tbut  in  a series  oi'  regions  of  length 
the  cutoff  plusiM  wavelength  X.  lmicI  electron  concentration  vaiiunce 
are  determined  from  some  onvlronmontnl  doscrii  tlon.  If  other  parameters 
such  ns  iixinl  electron  concentration  and  strlntlon  site  pariuiioter  are 
given,  the  formulations  of  Section  3.2  can  be  used  to  compute  o‘^  . dhesnut 
(and  TEMPO  following  his  work)  used  u profile  defined  ns  n^  oxp(-r*/o*)  , 

If  the  factors  accounting  for  the  difference  In  profile  doflnlllon,  the 
roLitionshlp  of  cutoff  siae  to  characteristic  siae  and  cutoff  wavolongth 
to  cutoff  size  arc  combined,  the  cutoff  wnvclengtli  should  be  taken  to  bo 
1.2  times  his  size  spoctnim  characteristic  parameter. 

5.1  Local  Effective  Strlat^on  Size 

The  local  effective  striution  size  can  be  taken  to  be 

0^,  - « O.Jib  . (19P) 

The  precision  of  this  relationship  should  be  20  to  SO  percent  as  estimated 
by  comparing  the  results  for  the  several  distributions  analyaed.  The  effec- 
tive size  is  defined  as  the  size  of  unifonii  strlations  In  un  environment 
having  the  same  electron  content  variance  needed  to  produce  the  same 
scattering.  The  same  effective  size  can  bo  used  for  relating  electron 
content  variance  and  scattering  angle  variance  (or  its  equivalent  phase 
derivative  variance)  and  for  relating  phase  variance  and  phase  derivative 
variance. 


5.2  Local  Scattering  Angle 

Several  directly  rolutod  parujneters  have  been  used  that  define 
the  hUbic  angular  scattering  or  vurlunco  of  transverse  derivative  of 
phase.  These  are  the  environmental,  plane  wave  scattering  angle  variance: 
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(200) 


^ I 


0 0 sin 

c 


a l2nj  • 


variance  of  transverse  derivative  of  phase: 


4Tt* 


°0  J 


and  variance  of  scattering  ns  viewed  by  radar: 


(201) 


(202) 


The  Inst  of  those  vaivinncos  indicates  the  proper  way  for  weighting  the 
scattering  as  a function  of  distance  from  the  radar, 


5.3 


by 


5.4 


Local  Phase  Variance 

The  locai  phase  variance  is  determined  from  the  phase  gradient 


^ of  “J.  • 

Integrated  Scattering  Angle 


(203) 


Inasmuch  as  any  of  the  local  scattering  angle  parameters  have  an 
integrated  value  that  is  a gaussian  variate,  the  variance  of  the  sum  is  the 
sum  of  tlio  variances  of  the  individual  contributions.  Thus, 


oj- 


I ‘ '1  ”i 


and 
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(204) 


(205) 
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5.5 


Integrated  Phase  Variance  and  Mean  Strlatlon  Size 


Tho  integrated  phase  also  becomes  u gausalan  variate  and  the 
variance  of  integrated  phase  will  bo 


0 


(206) 


One  would  like  to  be  able  to  obtain  the  integrated  phase  variance  from  the 
integrated  derivative  variance  by  use  of  a mean  offoctive  sire.  Thus, 


( 

1. 


(207) 


Comparison  of  equations  206  and  207  show  that  this  can  be  done  if  the 
mean  effective  size  is  formed  as 


i M 


i: 
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-M 


(208) 


(Because  they  are  related  by  a constant,  could  have  boon  used  instead 
of  In  liquation  208.)  If  the  mean  effective  strlution  size  is  deter- 
mined in  this  manner,  for  the  total  scattering  layer 
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?. 

'll 
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(209) 


5.6  Angular  Scattering  of  Power 


In  tho  limit  of  small  total  phase  perturbation  • there  is 

little  power  in  the  scattered  signal  and  its  angular  distribution  is  spread 
by  diffraction.  The  rms  scattering  angle  is 

tjj'  =«  / [1  ■ cxpC-o^)]  , (210) 

which  is  that  expected  for  diffraction  from  a structure  tlio  size  of  a 
correlation  length. 
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In  tho  limit  of  strong  scattering  (a^»l)  or  n uniform  strl- 
ation  weakly  scattering  plusmu,  the  power  is  gaussinnly  distributed  with  ii 
variance  given  by  liquation  210.  In  the  strong  scattering  multipath  case, 
the  number  density  of  paths  Is  gaussinnly  distributed  with  an  amplitude 
distribution  that  is  independent  of  the  distance  of  the  ray  from  the  direct 
line  of  sight.  Tor  weak  scattering  the  angular  power  spectrum  may  not  he 
gaussian  but  Its  variance  is  still  reasonably  approximated  by  lUiuntion  210 
and  usually  the  spectral  differences  are  only  significant  If  the  scattering 
is  Insignificant. 

5.7  Received  Power  Statistics 

If  the  signal  that  propagates  through  the  scattering  medium  l.s 
received  by  an  antenna,  the  received  power  (or  sum  channel  In  a mat)Dpul.so 
rocolvor.)  comprises  an  attenuated  but  unmodified  component  and  both  in- 
phase  and  quadrature  random  compone.its.  If  the  antenna  has  a sum  cluiniicl 
one-way  power  gain  (liaj  nn  angle  u from  the  beam  venter,  a .t-dll,  full 
width,  one-way  power  l)eaiiiwldth  of  , and  the  source  of  the  received 
energy  Is  an  angle  a from  the  hoaiii  center,  the  unmodified  component  has 
a power 

l'„  = I’o  e,\pl-o' ) btal  , (Jll  j 

where  I’j  Is  tlu;  power  that  would  be  received  In  the  absence  of  scattering 
If  tlic  antenna  has  been  pointed  at  tlie  source. 

If  a parameter  Is  defined  relating  the  scntterlng  to  the  atiteniui 

beamwldth , 

u 5,2.-  il'V'jJ  , (212J 
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the  rndur  components  onch  have  n voltaito  varinnce  (i.ei,  moan  power)  given 
hy 


^ [l-expc-oj)]  ^^exp 
“ ^ /ii+l  L b 


(213) 


If  the  antenna  Is  pointed  at  tho  source  so  that  a*0,  tho  power 
can  be  considered  to  bo  an  unmodified  signal  plus  a randomly  phased  noise 
component,  The  "signal-to-nolsc"  ratio  in  this  representation  Is 

S/N  y^l  c VU-e  ''’),  (214) 


It  can  bo  recognised  that  the  amplitude  distribution  will  be  Rician, 

5.8  Difference  Channel  Signal 


TItc  difference  channel  will  also  contain  both  an  unmodifiod  signal 
and  a modified,  scattered  component.  If  the  difference  channel  gain  is 
sln^ (p(^)  <Ua/2 ) , which  Is  normally  sero  fur  «»0,  the  unmodified  component 
iwhlch  is  in  iiuadrature  to  the  unmodified  sum  channel  component)  has  a 
liower 


= Pj  sln^Ipctl  C’lla/:) 


(21  SI 


Tho  total  power  In  the  modified,  random  component  will  be 
- I’o 


As 


ir\3/j  I'  -cxp(-«i)]  lUm'/d)  slnMpa 


whore  o'  « rtZ/T+'u/'T  and  p = l‘4/0|^  . 


' ) 1 + -tl 

ex' 


(21b) 
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5.9  L1in1t«t1ons 

The  analysis  has  assumed  that  sufficient  Independent  structures 
are  encountered  that  the  central  limit  theorem  results  in  gaussian  distri- 
butions of  integrated  phase  and  spatial  derivative  of  phase, and  the 
analysis  has  dropped  terms  that  uro  significant  only  in  the  antenna  near 
zone.  Both  these  assumptions  are  valid  for  most  problems  of  interest. 

The  angular  power  spectrum  is  not  necessarily  gaussian  if 
oj  « 1 but  ugnln  this  Is  not  usually  a regime  of  concern. 

Possible  correlations  between  the  sum  and  difference  channel 
signals  have  not  beon  e.\plorod.  Moreover,  the  specific  analysis  hero 
assumed  that  the  target  displacement  was  along  one  axis  of  the  monopulso 
difference  beam  and  the  striations  were  aligned  with  the  other  axis. 
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APPENDIX  A 
A USEFUL  INTEGRAL 


ConsUk'i’  ill!  t'ori 


ws 

I ■ J l’(x)  oxpf-{tix*  + bx  + c)]  ilx 


where  PCxI  is  u polynomial  in  x and  a,  b,  and  c nro  complex  constants 
(or  functions  of  variables  not  involved  in  the  Intcurnt  Ion ) , Tlu'  pi'ocosscs 
needed  to  perform  the  Integrations  are  straightforward  but  the  arithmetic 
is  often  tedious. 

The  procedure  is  to  use  a linear  transformation  of  variable  to 
remove  the  linear  term  in  the  exponent.  The  constant  term  can  then  be  taken 
outside  the  integration,  leaving  only  the  second  power  of  the  variable.  The 
polynomial  will  have  transformed  into  a new  polynoml.nl  of  the  same  order. 

The  integral  of  odd  order  terms  will  bo  zero  and  the  integral  of  oven  order 
terms  is  u well  known  definite  integral.  As  an  example,  lot  the  polynomial 
bo  a quadratic 

I’Cx)  » po  ♦Pi  X + ps  X*  (A2) 

Use  the  transformation 


* “ y ■ 2 «- 


This  produces 

I ■ / [(l’»  ‘ * i’"’ 

• 00 

X exp  ( -a  y»  ^ ^ - c)  ‘ly 
00 

/ b*  - 4ttC\  fL.  ..  >2L 

- oxp  ^ 5^  jJ{Po  2a 

•03 

This  final  Intcaval  has  the  value 


. ^ j y + pjy*j 
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